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ABSTRACT 


This  paper  considers  multiple  scattering  effects  of  millimeter 
wave  propagation  in  rain.  Using  the  multiple  scattering  theory  devel¬ 
oped  previously  by  the  authors  for  beam  and  plane  waves  propagating  in 
an  isotropically  scattering  discrete  random  medium,  rainfall  attenuations 
have  been  computed  for  intensities  ranging  from  0-100  mm/hr  at  wave¬ 
lengths  3-30  mm.  The  results  are  compared  with  existing  measurements  in 
the  literature.  The  good  agreement  between  theory  and  experiment  indi¬ 
cates  the  general  applicability  of  the  theory  to  rainfall  attenuation 
predictions  and  area  rainfall  intensity  estimations.  The  theory  also 
agrees  with  calculations  based  on  existing  theory  that  accounts  for  drop 
size  distributions.  The  difference  in  attenuation  of  plane  and  beam  waves 
propagating  in  rain  of  the  same  intensity  has  found  not  to  exceed  3  db/Km. 
The  limits  of  applicability  of  the  theory  have  also  been  indicated. 
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I.  INI AODUCTION 


In  recent  years,  there  has  been  increasing  interest  in  the  utili¬ 
sation  of  millimeter  wavus  for  communications  and  remote  sensing  [I.’]. 

This  frequency  spectrum  is  attractive  because  it  offers  the  promise  of 
large  bandwidth,  reliaf  to  the  already  crowded  radio  frequency  spectrum 
and  high  gain,  high  resolution  antennas  of  reduced  size  and  weight.  How¬ 
ever,  atmospheric  losses  impose  some  limitations  on  the  application  of 
millimeter  waves.  Millimeter  waves  are  attenuated  by  atmospheric  absorp¬ 
tion  and  scattering.  Except  for  the  1.35  cm  water  vapor  absorption  line 
and  the  S  mm  oxygen  absorption  line  [3  and  Fig.  1],  the  most  pronounced 
effect  is  due  to  rainfall. 

Attenuation  of  cm  and  non  waves  by  atmospheric  precipitations  was 
first  predicted  by  Stratton  [4],  in  a  theoretical  paper  published  in  1930, 

In  it  he  extended  an  earlier  investigation  of  Mie  [5],  and  concluded  that 
rainfall  might  attenuate  waves  of  5  cm  or  less  appreciably.  Since  then, 
there  has  been  numerous  work  done  both  experimentally  and  theoretically. 
However,  Medhurst  [6]  was  the  first  to  make  a  systematic  comparison  between 
the  then  available  experimental  data  in  the  literature  and  the  theoretically 
predicted  attenuation  versus  rain  rate.  In  the  paper,  Medhurst  extended 
the  work  of  Ryde  and  Ryde,  which  was  based  on  the  Mis  theory.  It  was  ob¬ 
served  that  in  many  cases ,  there  is  fairly  good  agreement  between  theory 
and  experiment,  but  in  some  cases  there  is  a  tendancy  for  the  measured 
attenuation  to  fall  well  above  the  maximum  theoretical  level.  In  the  years 
that  have  elapsed  since  Medhurst 's  paper,  several  experiments  have  been 
reported.  Again  they  exhibit  the  same  characteristics  as  Medhurst  ob¬ 
served. 

It  is  well  known  that  for  a  rain  of  given  rate  the  measured  attenua¬ 
tion  is  widely  scattered.  It  generally  varies  with  temperature,  humidity, 
wind  velocity  and  even  the  drop  size  distribution.  In  almost  all  of  the 
reported  measurements,  there  is  a  conspicuous  lack  of  records  of  the 
local  weather  conditions.  Because  of  this  inadequacy  in  meteorological 
information,  it  is  difficult  to  draw  any  definite  conclusion  as  to  the 
applicability  of  Mie  theory  to  rainfall.  However,  since  the  approach,  based 
on  Mie  theory,  merely  considers  the  interaction  of  a  plane  wave  with  a 
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single  isolated  scatterer,  and  the  total  attenuation  is  simply  to  be 
the  additive  sum  of  absorption  and  scattering  from  each  individual 
scatterer,  it  is  possible  that  this  single  scattering  theory  is  not  suf¬ 
ficient  to  describe  an  wave  propagation  in  rain.  In  fact.  Crane  [7]  has 
recently  shown  that  it  i3  important  to  consider  the  multiple  scattering 
effects  of  rain  for  frequencies  above  20  GHz. 

The  formulation  of  multiple  scattering  effects  of  a  general  electro¬ 
magnetic  wave  propagation  in  an  isotropically  scattering  random  medium 
was  given  by  the  authors  [81.  These  results  give  explicit  expressions 
for  the  beam,  spherical  and  plane  sources  as  function  of  density,  cross 
sections  and  source  parameters.  It  is  shown  that  the  beam  wave  average 
intensity  is  essentially  coherent  for  most  practical  applications,  whereas 
the  plane  wave  case  has  significant  contribution  from  both  coherent  and 
incoherent  propagation.  In  general,  the  effect  of  multiple  scattering  is 
to  increase  the  received  power,  therefore  the  predicted  attenuation  is 
less  than  that  of  the  single  scattering  theory. 

The  purpose  of  this  paper  is  to  apply  these  results  to  mm  wave  prop¬ 
agation  in  rain.  It  is  our  intention  to  compare  with  the  existing  theory 
and  reported  measurements.  The  good  agreement  with  published  experimental 
measurements  indicates  the  applicability  of  the  theory  to  rainfall  attenua¬ 
tion  predictions.  The  results  are  of  interest  to  the  electrical  engineers 
concerned  with  the  utilization  of  millimeter  wave  bands  for  radio  relay, 
terrestrial  and  satellite  comnunication  systems.  To  the  meteorologist, 
the  results  are  useful  for  estimating  area  rain  intensity  from  attenuation 
measurements. 


II.  THEORY  Of  SCATTERING  PROPAGATION 

In  this  section,  ue  summarize  some  results  obtained  previously  [0], 
These  results  will  be  used  to  predict  rainfall  attenuation  of  mm  wave 
propagation. 

A.  SINGLE  DROP 

We  consider  a  plane  wave  (say  component  x,  of  the  electric  fiald) 


(1) 
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incident  on  a  rain  drop  at  r  s  C,  where  k  =  2m/X  is  the  free  space  wave 
number  and  e'*wt  time  dependence  is  assumed.  The  scattered  wave  in  the 
far  field  is 


f  (?,£> 


where  f  is  the  scattering  amplitude,  and  z  and  6  are  directions  of  inci¬ 
dence  and  observation  respectively.  From  the  optical  theorem  we  have 


0*  *<M> 


o  =  /  |f<s.*>rdfi  (4) 

*  4ir 

where  o  and  o<  are  total  and  scattering  cross  sections  respectively. 

B„  HIE  THEORY 

A  complete  theory  of  scattering  of  a  plane  wave  by  a  sphere  of  any 
composition  in  an  external  medium  was  obtained  by  Hie  [5]  in  1908. 

This  has  become  the  oasis  for  most  subsequent  discussions  in  the  sub¬ 
ject.  Hie  separated  the  field  equations  in  spherical  coordinates  into 
two  groups  -  one  transverse  magnetic  and  the  other  transverse  electric. 
*!atchirg  the  boundary  conditions  at  the  surface  of  the  sphere,  he  obtained 
the  cross  sections  of  a  spherical  particle. 


(-ke)  I  (2n+i)(a  +  b  ) 

,  n  n 

n=l 


l  (2ntl)(|a  |2  ♦  |bj2> 


where  and  are  spherical  Bessel  functions  of  orde:  n,  and  they  are 
given  in  terms  of  the  complex  index  of  refraction,  the  radius  of  the 
sphere  a,  and  wavelength  of  the  medium  X. 

The  total  cross  section  o  and  scattering  cross  section  o  given  by 
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where  is  the  positive  root  of 

„  tanh-1a 
a  o 


and 


g(u,  o Jo)  = 
8 


°m  .1  0  TO«U  0 

u  tanh  u)  +  (-^|-) 


(1  - 


(10) 


(11) 


Eq.  (9)  is  obtained  through  a  Fourier  transform  technique  applied  to 
the  integral  equation. 

-do|r-f'| 

<T(r)>  =  C(r)  ♦  d  /  - - =-  if(f',2,)r  <T(r')>  dV(r')  (12) 

if-f’i2 

It  Is  noted  that  as  the  wave  propagates  in  the  medium,  the  coherent  in¬ 
tensity  decreases  while  the  incoherent  intensity  increases  in  relative 
magnitude.  In  addition,  as  described  elsewhere  [8],  the  average  intensity 
of  a  beam  wave  propagating  in  a  discrete  random  medium  is  essentially 
coherent,  and  it  is  approximately  given  by 


<T  >  «* 
3e  ^ 


exp 


1  t 


doz 


|K| 


2  2 
“Sr 

o 


w 


2  4 

W 

o 


-doz  - 


2W 


2  + 


doz 

k|2w2 
'  o 


1+ 


doz 

K  |2vP 


1  o 


(13) 


Therefore,  beam  wave  experiences  heavier  attenuation  as  it  propagates 
through  rain  than  plane  wave.  However,  the  high  gain  associated  with 
a  narrow  beam  in  practice  will  more  than  compensate  for  this  additional 
loss. 


III.  APPLICATIONS  TO  RAINFALL  ATTENUATION 
A.  INDEX  Of  REFRACTION 

The  refractive  index  m  *  n  +  ix  of  water  is  generally  complex  at 
ram  wavelengths  and  it  is  a  function  of  frequency  and  temperature. 
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Saxton  111]  has  calculated  curves  showing  the  variation  of  the  refractive 
index  of  water  over  the  wavelength  range  2  m.~  to  10  cm,  which  covers 
practically  the  whole  dispersion  region.  He  also  described  the  effect  of 
temperature.  In  addition,  the  calculated  values  were  compared  with  ex¬ 
perimental  data  obtained  by  other  independent  investigators.  They  were 
found  to  agree  very  well.  In  Table  I,  we  give  the  refractive  index  of 
water  for  the  wavelengths  under  investigation.  Values  at  20°C  are  taken 
from  Saxton,  and  those  at  18°C  are  from  Lowan  [9],  The  values  at  10°C 
and  0°C  are  those  given  by  Herman,  Browning  and  Battan  [10]. 

B.  DROP  SIZE  DISTRIBUTION  AND  TERMINAL  VELOCITY 

In  order  to  relate  accurately  the  attenuation  to  an  actual  rain  of 
given  intensity,  it  is  necessary  to  know  the  droo  size  distribution. 
However,  natural  rainfall  varies  with  prevailing  wind,  temperature,  and 
other  conditions.  Furthermore,  for  a  given  propagation  path,  the  rainfall 
is  seldom  uniform.  It  would  seem  practical  then  to  compute  the  rainfall 
attenuation  for  rain  with  uniform  average  droj  size.  For  our  computations 
we  choose  the  median  drop  size  from  the  Laws  and  Parsons  distribution  [12]. 
The  median  drop  divides  the  larger  and  smaller  drops  of  a  given  rair.  in¬ 
tensity  into  two  groups  of  equal  volume.  This  median  dron  radius  is  s'pown 
in  Fig.  2,  for  rainfall  intensities  ranging  from  0  to  100  mrn/hr. 

The  terminal  velocities  of  water  drops  in  still  air  are  given  bv 
Gunn  and  Kinzer  at  a  pressure  of  760  mm  and  temDerature  of  20°C,  and  ere 
readily  available  [13].  This  is  shown  in  Fig.  3. 

It  should  be  mentioned  that  in  actual  rainfall  there  is  usually  a 
vertical  component  of  wind  associated  with  it.  Even  though  it  mav  be 
small,  the  vertical  component  is  very  important.  It  is  the  up-and-down 
motion  of  air  particles  that  is  principally  responsible  for  much  of  the 
atmospheric  activities. 

C .  THEORETICAL  ATTENUATION 

For'  comparison  with  some  of  the  readilv  available  measurements  in 
the  literature,  using  Eas.  (9)  and  (13),  numerical  computations  have  been 
made  for  attenuation  at  wavelengths  3,  5,  9,  12.5,  17,  20  and  30  mm  for 
rainfall  intensities  ranging  from  0  to  100  mm/hr,  using  a  CDC-64G0  com¬ 
puter.  According  tc  Humphrey's  classification  of  rain  [14],  the  inten- 
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TABLE  I 

Variation  of  Refractive  Index  of  Watar  ro  =  n  +  ix 
with  Wavelength  and  Te.oeraturg 


X(n»n) 


3 

5 

9 

12.5 

17 

20 

30 

2.759 

3.192 

4.681 

5.837 

6.929 

3.106 

3.767 

5.740 

6.943 

7.851 

3.410 

4.210 

5.550 

6.410 

7.200 

7.560 

8.180 

3.509 

4.370 

5.7B0 

6.530 

7.230 

7.643 

8.251 

1.241 

1.766 

2.732 

3.005 

2.929 

1.663 

2.219 

2.920 

2.808 

2.397 

1.940 

2.510 

2.850 

2.860 

2.650 

2.550 

1.960 

2,010 

2  524 

2.850 

2.770 

2.570 

2.395 

1.827 
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sities  covered  encompass  rates  of  precipitation  up  to  the  grades  of  ex¬ 
cessive  rain  and  cloudburst  encountered  in  nature. 

The  total  and  scattering  cross  sections  are  the  accurate  values 
computed  by  Lowan  using  Hie  theory.  Cross  sections  at  X  =  12.5  and 
20  mm  are  interpolated  values  from  Lowan' s  data.  Curves  for  these  cross 
sections  versus  water  drop  radius  at  the  wavelengths  we  mentioned  are 
shown  in  Figs.  (4a)  and  (4b). 

In  practice,  it  has  been  customary  and  convenient  to  express  rain 
attenuation  as  a  function  of  rainfall  rate,  p  in  mm/hour.  This  factor 
is  proportional  to  the  terminal  velocity,  which  in  turn  depends  on  the 
size  of  drops.  In  order  to  get  a  quantitative  measure  of  the  attenuation 
as  a  function  of  rainfall  rate,  we  must  have  a  relationship  between  p  and 
"the  average  number  of  drops  per  cubic  meter  d  (density  of  distribution). 
This  is  given  by 

p  =  15.1  vda3 

where  v  is  the  terminal  velocity  in  meters  per  second  and  a  is  the  rain¬ 
drop  radius  in  centimeters. 

(Jaing  the  values  of  a,  v,  o  and  o  shown  in  Figs.  2-4,  theoretical 

s 

attenuation  curves  are  calculated  and  plotted  in  Figs.  5-11.  In  addition, 
in  each  figure  is  shown  the  Hie  theory  curve  based  on  Laws  and  Parsons 
drop  size  distribution,  along  with  the  Hie  xheorv  curve  of  a  hypothetical 
rain  of  uniform  drop  size  specified  by  the  median  of  Laws  and  Parsons 
distribution.  The  Hie  theory  cases  are  taken  from  Hedhurst. 

D.  HFASURED  ATTENUATION 

The  attenuation  of  electromagnetic  radiation  by  rainfall  may  be 
measured  in  two  ways.  In  the  first  method,  transmitting  and  receiving 
antennas  are  set  up  at  the  opposite  ends  of  a  line-of-sight  path.  In  the 
second,  reflecting  targets  are  placed  on  the  line  of  radar  sight.  Echo 
returned  from  the  reflector  is  used  to  measure  the  attenuation.  In  both 
cases  the  rain  intensities  are  estimated  from  rain  gauges  located  along 
the  path.  There  is  however  no  uniformity  of  rain  gauges  used  for  the 
rain  intensities.  Although  a  great  majority  of  the  experiments  referred 
to  in  this  paper  used  funnel  and  gradual  type  rain  gauges,  they 
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are  not  all  of  equal  quality.  Some  are  quits  advanced  devices,  others 
are  rather  primitive.  Hence,  it  is  difficult  to  compare  measurements 
ty  different  workers.  Furthermore,  the  degrees  of  reliability  vary 
according  to  the  number  of  rain  gauges  and  path  lengths. 

An  apparent  and  serious  failure  of  almost  all  previously  published 
experiments  is  the  neglect  of  local  climate  conditions.  Seldom  Is  the 
ambient  temperature,  prevailing  wind  and  humidity  recorded.  Considera¬ 
tions  of  these  meteorological  parameters  are  important,  because  raindrop 
absorption  and  scattering  depend  c.i  the  complex  refractive  indu  ;s,  which 
are  functions  of  temperature.  Since  *  vertical  component  of  wind  may  in¬ 
crease  or  decrease  the  terminal  velc>-  .y  of  raindrops,  o  condition  can 
be  set  up  where  the  density  of  raindrops  in  the  propagation  path  may  be 
different  from  what  is  indicated  by  the  rain  gauges.  Furthermore,  the 
presence  of  wind  may  complicate  the  situation  by  giving  a  spatial  distribu¬ 
tion  of  drops  such  that  smaller  drops  are  carried  further  downwind  than 
large  drops.  Lastly,  the  presence  of  atmospheric  gasses  and  water  vapor 
may  modify  the  attenuation  rate  in  some  instances.  Table  II  shows  a 
summary  of  the  parameters  of  rainfall  attenuation  measurements  in  tne 
literature  considered  in  this  paper.  [15j-[19],  [21]-[2*0,  [27] 

For  comparison  with  the  theoretically  predicted  attenuations,  the 
experimentally  measured  values  are  shown  in  Figs.  5-11,  along  with  the 
theoretical  curves. 

E.  DISCUSSION 

From  the  graphical  presentation,  it  is  seen  that  the  multiple  scat¬ 
tering  prediction  of  attenuation  of  mm  wave  propagation  in  rain  of  uni¬ 
form  drop  size  specified  by  the  median  is  in  good  agreement  with  experi¬ 
mental  measurements  and  Mie  theory,  when  the  latter  assumes  a  drop  size 
distribution  given  by  Laws  and  Parsons,  down  to  5  mm  range.  For  the 
cases  of  3  and  5  mm  shown  in  Fig.  5  and  6,  the  theoretical  values  deviate 
considerably  from  both  experimental  and  Mie  for  higher  rainfall  rates. 

This  is  more  pronounced  in  the  3  mm  case.  However,  at  low  intensities 
the  attenuations  are  the  same  for  theory  and  experiment.  These  are  not 
inconsistent  with  the  present  theory.  For,  at  these  wavelengths,  the  drop 
sizes  are  increasingly  comparable  to  wavelength  as  the  rain  intensities 
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get  higher,  the  isotropic  scattering  approximation  is  no  longer  justified 
and  the  theory  does  not  apply,  whereas  for  lover  rain  intensities  drop 
sizes  are  still  small  compared  to  wavelengths. 

The  situation  at  12.5  mm  displayed  in  Fig.  8  is  of  a  different  char¬ 
acter.  The  only  available  experimental  data  (the  reason  may  be  12.5  mm 
is  generally  considered  to  be  an  atmospheric  stop-band  for  mm  wave)  fall 
above  the  theoretically  predicted  attenuation  by  a  considerable  amount. 

The  discrepancy,  however,  may  not  be  due  to  the  inapplicability  of  the 
theory  but  rather  due  to  the  inadequacy  in  the  handling  of  meteorological 
information.  Aside  from  the  fact  that  additional  attenuation  is  experi¬ 
enced  from  the  proximity  of  12.5  mm  wavelength  to  the  1.35  nan  water  vapor 
absorption  line  (see  Fig.  1),  there  are  evidences  [25]  that  Hawaiian  oro- 
graphical  rains  generally  have  different  structures,  therefore  the  rain¬ 
drop  size  distribution  is  not  the  same  as  that  measured  by  Laws  and  Parsons 
which  was  used  for  the  theoretical  computation.  Furthermore,  as  Anderson, 
Day,  Fr?res  and  Stokes  [20]  have  observed,  there  was  a  light  wind  associ¬ 
ated  with  the  orographic  precipitation.  However,  this  was  not  taken  care 
of  in  their  interpretation  of  the  experiment. 

It  is  apparent  from  this  study  that  attenuation  does  not  increase 
linearly  with  rainfall  rate.  The  only  possible  exception  of  this  is  around 
9  mm. 

In  light  of  the  high  degree  of  collimation  obtainable  from  mm  wave 
beam  u*ing  moderate  size  antennas,  it  is  of  interest  to  investigate  the 
effects  of  wave  types.  In  particular,  we  have  studied  plane  and  beam  waves. 
As  shown  by  Eqs.  (9),  (13)  and  elsewhere  [8],  the  beam  wave  average  inten¬ 
sity  is  essentially  coherent.  Furthermore,  the  coherent  intensity  normal¬ 
ized  to  free  space  is  the  same  as  plane  wave.  This  has  been  computed  and 
shown  in  Figs.  5-11.  On  the  other  hand,  plane  wave  has  significant  con¬ 
tribution  from  both  coherent  and  incoherent  intensities  as  shown  in  Figs. 
5-11.  Within  the  range  of  applicability  of  the  present  multiple  scattering 
theory,  the  maximum  difference  between  the  two  wave  tvpes  is  about  3  db/Km. 

Physically  this  corresponds  to  the  fact,  because  of  the  finite  beam 
size,  waves  scattered  out  of  the  beam  which  would  normally  be  multiply 
scattered  and  eventually  reach  the  receiving  antenna,  as  in  the  plane  wave, 
will  not  do  so  for  beam  wave  until  further  away.  By  then,  the  coherent 
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intensity  is  attenuated  so  much  that  they  both  appear  as  r.nise  to  the 
receiver. 

The  effects  of  non-sphericity  to  raindrops  have  not  been  considered 
in  this  paper.  It  has  been  shown  that  large  falling  raindroDS  tend  to  be 
flattened  at  the  bottom.  Therefore,  it  is  plausible  that  rainfall  atten¬ 
uations  are  polarization  dependent.  Semplak  [26]  has  recently  shown  that 
there  is  a  significant  difference  in  mm  wave  attenuations  for  horizontal 
and  vertical  polarizations  for  the  same  rainfall.  It  is  also  clear  from 
the  measurements  cited  above  that  there  are  many  .mbigulties  concerning 
rain  attenuation  experiments.  It  would  seem  desirable  and  important  to 
carry  out  a  series  of  experiments  with  due  care  exercised  in  the  measuring 
and  recording  of  all  the  pertinent  meteorological  Darameters.  The  addition¬ 
al  experiments  will  go  a  long  way  to  elucidate  the  Dresent  uncertainties. 


IV.  CONCLUSIONS 

This  work  shows  the  effect  of  multiple  scattering  on  mm  wave  cropaga- 
tion  in  rair..  It  has  been  shown  that  using  the  median  dron  size  given 
bv  Laws  and  Parsons  distribution,  computed  theoretical  attenuation  agree 
well  with  measurements  at  wavelengths  range  from  S  mm  and  un  for  all 
grades  of  rainfall  intensities.  Furthermore,  the  multiple  scattering 
theory  is  in  agreement  with  the  Mie  theory,  when  the  latter  assumes  a 
full  spectrum  of  drop  sizes.  In  addition,  the  effect  of  wave  type  have 
been  investigated  to  show  a  maximum  difference  of  3  db/Kro  for  nlane  and 
beam  waves.  The  limits  of  applicability  of  the  theory  have  also  beer,  in¬ 
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